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Abstract 

It is well known that three momentum wheel actuators can be used to control the attitude 
of a rigid spacecraft and that arbitrary reorientation maneuvers of the spacecraft can be 
accomplished using smooth feedback. If failure of one of the momentum wheel actuators 
occurs, we demonstrate that two momentum wheel actuators can be used to control the atti- 
tude of a rigid spacecraft and that arbitrary reorientation maneuvers of the spacecraft can be 
accomplished. Although the complete spacecraft equations are not controllable, the spacecraft 
equations arc small time locally controllable in a reduced nonlinear sense. The reduced 
spacecraft dynamics cannot be asymptotically stabilized to any equilibrium attitude using a 
time-invariant continuous feedback control law, but discontinuous feedback control strategies 
are constructed which stabilize any equilibrium attitude of the spacecraft in finite time. Conse- 
quently, reorientation of the spacecraft can be accomplished using discontinuous feedback 
control. 


* Please send all correspondence to Professor N. Hams McClamroch, Department of 
Aerospace Engineering, University of Michigan, Ann Arbor, MI 48109-2140. 
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1. Introduction 

We consider the attitude control of a spacecraft modeled as a rigid body. It is well 
known that three actuators, either gas jets or momentum wheels, can be used to control the 
attitude of a rigid spacecraft and that arbitrary reorientation maneuvers of the spacecraft can 
be accomplished using smooth feedback 1-7 . If failure of one of the actuators occurs, then one 
is left with only two actuators. In this paper, the attitude stabilization problem of a rigid 
spacecraft using only two control torques supplied by momentum wheel actuators is con- 
sidered. Since we are considering a space-based system, the problem considered here, namely, 
the attitude stabilization of a spacecraft operating in an actuator failure mode, is an important 
control problem. It is assumed that the center of mass of the system consisting of the space- 
craft and the momentum wheel actuators is fixed in space. 

Attitude stabilization of a rigid spacecraft using two momentum wheel actuators is not a 
mature subject in the literature. Controllability results for a rigid spacecraft controlled by 
momentum wheel actuators are presented in Ref. 8. We mention that most of the previous 
researchers have considered the problem of controlling a rigid spacecraft using two gas jet 
actuators 8-22 . Attitude stabilization of a rigid spacecraft using two gas jet actuators is con- 
sidered in Refs. 8-13. Refs. 14-22 consider only the stabilization of the angular velocity 
equations of a rigid spacecraft using two gas jet actuators. 

We consider the attitude stabilization of a spacecraft using control torques supplied by 
two momentum wheel actuators about axes spanning a two dimensional plane orthogonal to a 
principal axis of the spacecraft. The linearization of the complete spacecraft dynamic equa- 
tions at any equilibrium attitude has an uncontrollable eigenvalue at the origin. Consequently, 
controllability and stabilizability properties of the spacecraft cannot be inferred using classical 
linearization ideas. The complete spacecraft dynamics is, in fact, not controllable. Under a 
rather weak assumption, the spacecraft dynamics is small time locally controllable at any 
equilibrium attitude in a reduced nonlinear sense. The reduced spacecraft dynamics cannot be 
asymptotically stabilized to any equilibrium attitude using time-invariant continuous feedback. 
Nevertheless, two different discontinuous feedback control strategies are constructed which 
achieves reorientation of the spacecraft in finite time. Using the concept of geometric phase 23 , 
a discontinuous feedback control strategy is presented based on the nonholonomic control 
theory in Ref. 24. An alternate discontinuous feedback control strategy, based on the fact that 
rigid body rotations do not commute, is also presented. 
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This paper is based on our ealier work presented in Ref. 10 and is a companion to Ref. 
11 and Ref. 12, which treat the attitude stabilization of a rigid spacecraft using two gas jet 
actuators. 


2. Kinematic and Dynamic Equations 

The orientation of a rigid spacecraft can be specified using various parametrizations of 
the special orthogonal group 5(9(3). Here we use the Z-Y-X Euler angle convention for 
parametrizing the orientation of the rigid spacecraft 25 . The corresponding rotation matrix is 
denoted as R (y, 0/ <J>), where y, 0, 4> are the Euler angles. We assume that the Euler angles 
are limited to the ranges -n < y < n, -n/2 < 0 < jc/2, -jc < <j> < k. Suppose tOj, to* CO3 are 
the principal axis components of the absolute angular velocity vector co of the spacecraft 
Then we have 25 

<j> = <»! + (02sin<|> tan0 + o^coscj) tan0, (2.1) 

0 = ©2COS<J> - ( 03 sin<}>, (2.2) 

y = 0)2sin<j) sec0 + cojcos^ sec0. (2.3) 


Next we consider the dynamic equations which describe the evolution of the angular 
velocity components of the spacecraft. Consider two momentum wheel actuators spinning 
about axes defined by unit vectors b x , b 2 fixed in the spacecraft such that the center of mass 
of the i -th wheel lies on the axis defined by b t , and a control torque - u, is supplied to the 
i-th wheel about the axis defined by b t by a motor fixed in the spacecraft Consequently, an 
equal and opposite torque u t is exerted by the wheel on the spacecraft We assume that 6, 
defines a principal axis for the i'-th wheel which is symmetric about b x . Further b x and b 2 
span a two dimensional plane which is orthogonal to a principal axis of the spacecraft and, 
without loss of generality, b t are assumed to be of the form 


bi = (b n , b i2 , Of, i = 1, 2. 


(2.4) 


The mass of spacecraft, wheel 1 and wheel 2 are denoted as m lt m 2 and m 3 respectively, 
and p lt p 2 , p 3 denote the position vectors of the center of mass of the spacecraft wheel 1 and 
wheel 2 respectively with respect to the center of mass of the whole system. Thus from the 
location of the wheels 


p 2 - Pi + d x b 1, 


(2.5) 
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P3 — Pi ^2^2« (2.6) 

where d\ t d 2 are constants. Since, by the definition of center of mass, 

3 

X"*;Pi=0, (2.7) 

«=l 

further manipulation of equations (2.5)-(2.7) gives expressions for p lf p 2 and p 3 which we 
denote as p,- = (c, i, c i2 , 0 ) T , i = 1, 2, 3. The total angular momentum vector of the system is 
given, in the spacecraft body frame, by 

R (y, 0, ty)H = /co + v , (2.8) 

where 


J = 


+ hi + 

3 

I (A - 



(2.9) 



1=1 

«= 2 






cn 

~ c il c i2 

0 



h - 

= m t 

~<i l c x2 

r 2 
c i 1 

0 

,i=l, 2, 3, 

(2.10) 



0 

0 

r 2, r 2 
W 1 ' 2. 



h 

= bib\ji. 




(2.11) 

h 

= b 2 b\j 2 . 




(2.12) 

V = 

- /^(co + b i®}) 

3 

+ 

+ b 2 Q 2 ), 


(2.13) 


where I / 2 , and / 3 denote the inertia tensors of the spacecraft, wheel 1 and wheel 2 respec- 
tively, j l is the moment of inertia of wheel 1 about the axis defined by b lt j 2 is the moment 
of inertia of wheel 2 about the axis defined by b 2 , and 0 lt 0 2 are the angles of rotation of 
wheel 1 and wheel 2 about the axes defined by b x and b 2 respectively. Here H denotes the 
angular momentum vector of the system expressed in the inertial coordinate frame. The angu- 
lar momentum vector H is a constant since there is no external moment about the center of 
mass of the system. Suppose Uj and u 2 are the control torques; then 

v = - (biUi + b 2 u 2 ). (2.14) 

Differentiating (2.8) with respect to time we obtain 

7d) = S (co)R (\y, 0, (j ))H + b^i + b 2 u 2 , (2.15) 


where 
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S(co) = 


- 

0 

-0)3 

0)2 


0)3 

0 

-0)j 


- 0)2 

0), 

0 


Note that 


1 1 - diag(/ |j, / 12» ^ 13 ). 

/ 2 = block diag(/ 2 j, / 22)* 

/ 3 = block diag(/ 31 , / 32 ), 

where /21, /31 are invertible 2x2 matrices, / n, / 12* f i3> f 22 ’ ^ 32 are nonzero real numbers 
and therefore / is a positive definite matrix of the form 

J = block diag (J ^2). ( 2 - 16 ) 

where J\ is an invertible 2x2 matrix and J 2 is a nonzero real number. 


3. Controllability and Stabilizability Properties 
In this section we consider the controllability and stabilizablity properties of the space- 
craft dynamics controlled by two momentum wheel actuators. Define 
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From Section 2 the complete spacecraft dynamics can be rewritten as 


V 

0 (lx2) 


0(2x1) 


S(G))/?(y, 9, <)>)// 


+ 


“1 

«2 • 
0 


( 3 . 1 ) 


<j) = CDj + o^sin^ tan 0 + CO3COS <J> tan 0 , ( 3 - 2 ) 

0 = oncost)) - 0)3sin<J>, ( 3 - 3 ) 

y = cosine}) see© + a^cosij) sec0, ( 3 - 4 ) 

where H is a constant vector. 

The linearization of the complete spacecraft dynamic equations ( 3 . 1 )-( 3 . 4 ) at any equili- 
brium attitude has an uncontrollable eigenvalue at the origin. Consequently, the controllability 
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and stabilizability properties of the complete spacecraft dynamics cannot be inferred using 
classical linearization ideas. However, from equations (2.4), (2.11)-(2.13) and the definition 

c = (0, 0, 1) T , (3.5) 

we have c T v = 0. Therefore from equation (2.8) we have 

c r /?(y, 0, <>)// = c T J to. (3.6) 

Since H is a constant vector, this equation represents a constraint on the motion of the space- 
craft irrespective of the controls applied. Thus the complete spacecraft dynamics is not com- 
pletely controllable'. Therefore we ask the following question: what restricted control and sta- 
bilization properties of the spacecraft can be demonstrated in this case? Our analysis begins 
by demonstrating that, under an appropriate restriction of interest, the spacecraft equations 
have restricted controllability and stabilizability properties. 

Consider equations (3.1)-(3.4) and suppose the angular momentum vector H of the sys- 
tem is zero. From equations (2.16), (3.5) and (3.6) it follows that the angular velocity com- 
ponent of the spacecraft about the uncontrolled principal axis is identically zero, i.e., co 3 a 0. 
Under such a restriction, the reduced spacecraft dynamics are described by 


<*>! = «!, 

(3.7) 

= w 2> 

(3.8) 

<J> = «! + G^sintJ) tan0. 

(3.9) 

0 = a>2Cos<J>, 

(3.10) 

V = co 2 sin<J> sec0. 

(3.11) 


Notice that the linearization of the equations (3.7)-(3.11) at any equilibrium has an uncontroll- 
able eigenvalue at the origin. Therefore analysis of the controllability and stabilizability pro- 
perties of the reduced spacecraft dynamics requires inherently nonlinear techniques. The fol- 
lowing results follow directly based on an analysis similar to that in Ref. 24. 

Theorem 3.1: The reduced dynamics of a spacecraft controlled by two momentum wheel 
actuators as described by equations (3.7)-(3. 11) are small time locally controllable at any 
equilibrium. 

Theorem 3.2: The reduced dynamics of a spacecraft controlled by two momentum wheel 
actuators as described by equations (3.7)-(3.11) cannot be asymptotically stabilized to any 
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equilibrium using a rime-invariant continuous feedback control law, but the reduced dynamics 
can be asymptotically stabilized to any equilibrium using a piecewise continuous feedback 
control law. 

Theorem 3. 1 follows from the fact that a sufficient condition for small time local control- 
lability given in Ref. 26 is satisfied by the equations (3.7)-(3.11). The first part of Theorem 
3.2 follows from the fact that a necessary condition for the existence of a time-invariant con- 
tinuous feedback control law given in Ref. 17 is not satisfied by equations (3.7)-(3.11); the 
second part is a consequence of small time local controllability 26 . The implications of the 
properties stated above are as follows. Suppose the angular momentum vector H is zero. 
Then the spacecraft controlled by two momentum wheel actuators can be controlled to any 
equilibrium attitude but the feedback control law must necessarily be discontinuous. Thus 
arbitrary reorientation of the spacecraft can be achieved under the restriction H = 0; If H * 0, 
equation (3.6) implies that reorientation of the spacecraft to an equilibrium attitude cannot be 
achieved. 


4. Feedback Stabilization Algorithms 

We restrict our study to the class of discontinuous feedback controllers in order to asymptoti- 
cally stabilize the reduced spacecraft dynamics described by state equations (3.7)-(3.11). 
Clearly, traditional nonlinear control design methods are of no use since there is no general 
procedure for the design of a discontinuous feedback control. However, an algorithm generat- 
ing a discontinuous feedback control which asymptotically stabilizes an equilibrium can be 
constructed, as suggested by the controllability properties of the system. Without loss of gen- 
erality, we assume that the equilibrium to be stabilized is the origin. We present two different 
discontinuous control strategies which stabilize the origin of equations (3.7)-(3.11) in finite 
time. 

4.1. Feedback stabilization based on nonholonomic control theory 
Consider a diffeomorphism defined by 
y ! = cos<|> /n(sec0 + tanQ) + ysin<}>, 
y 2 = oc^sece -y^s, 

y 3 = <J>> 


(4.1) 

(4.2) 

(4.3) 
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y 4 = coj + C 02 sin<|) tan 0 , 

( 4 . 4 ) 

y 5 = sin<J> ln(sec 0 + tan 0 ) - \j/cos<j>, 

( 4 . 5 ) 


If we now define the feedback relations 


“2 


— sin<J>sin 0 (1 - y 5 sin<|>sin0) 
COS0 y 5 COS0 


v i 

v 2 


costyy ssin^sec 2 ©^ 
0 


— yfy i + cos<J)(sec 0 tan 0 o )2 - y 5 y 4 tan 0 <» 2 ) 
cos 4 >(y 4 tan 0 ci >2 + sin^sec^co^) 


( 4 . 6 ) 


then the reduced spacecraft dynamics (3.7)-(3.11) are described in the new variables by the 
normal form equations 


II 

K> 

( 4 . 7 ) 

y 2 = v i. 

( 4 . 8 ) 

y 3 = > 4 . 

( 4 . 9 ) 

Vi. 

II 

< 

K) 

( 4 . 10 ) 

ys-y^v 

( 4 . 11 ) 


From equations (4.1)-(4.5), notice that co 1 = to 2 = <l> = 0 = ¥ :=0 implies that y l = y 2 = y 3 = 
y 4 = y 5 = 0. Hence asymptotic stabilization of equations (3.7)-(3.11) to the origin is 
equivalent to asymptotic stabilization of the normal form equations (4.7)-(4.11) to the origin; 
hence we consider asymptotic stabilization of the normal form equations. The normal form 
equations (4.7)-(4.11) are in a familiar form which has been studied in Ref. 24 and therefore 
can be stabilized by the following discontinuous control strategy. 

• First, transfer the initial state of the normal form equations (4.7)-(4.11) to the equilibrium 
state ( 0 , 0 , 0 , 0 , y$ ), for some y$ , in finite time. 

• Next, traverse a closed path y in the (y j, y 3 ) space in finite time, where the path y is 
selected to satisfy 

-y} =J t 3 , i^3; 

this transfers the state ( 0 , 0 , 0 , 0 , y 5 ) to the origin in finite time. 


( 4 . 12 ) 
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Here we consider a rectangular path y in the (y lt y 3 ) space formed by line segments from 
( 0 , 0 ) to (y i , 0 ), from (yj, 0 ) to (y\,y]), from (y\,y* 3 ) to ( 0 ,y^), and from ( 0 ,y^) to 
(0, 0). For such a path, the line integral in equation (4.12) can be explicidy evaluated as y[y^ 
so that equation (4.12) becomes 

-ys 1 =ylyl’ (413) 

and the parameters y* and y 3 specifying the particular rectangular path are chosen to satisfy 
the above equation. 

Throughout, assume k > 0, and define 

x 2 lx 2 l 


G(x lt x 2 ) = 


k if {*! + 
Ui + 

-k if {xj + 


2k 

x 2 \x 2 \ 

2k 

x 2 \x 2 \ 


2k 


> 0 } or 

= 0 and x 2 > 0 } 
- < 0 } or 


x 2 \x 2 \ 

{x, + — — — =0 and x 2 < 0 ) 
2Jc 


0 if (X[ = 0 and x 2 = 0 } 

We use the well-known property that any initial state of the system 


x 1 = x 2 , 

x 2 = - G(xj -Xj, x 2 ), 

is transferred to the final state (T[, 0 ) in a finite time. 

We now present a specific feedback control algorithm which stabilizes the spacecraft to 
the origin in finite time; this feedback control algorithm implements the approach just 
described. 


Maneuver 1: Apply 

vi = - G(y lf y 2 ) , 
v 2 = -G(y 3 ,y 4 ) , 

until (y j, y 2 , y 3 , y 4 , y 5 ) = (0, 0, 0, 0, y} ) where y 5 ‘ is arbitrary; then go to Maneuver 2. 
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Maneuver 2 : If >5 £ 0, choose y\ = -y 3 = VO’S 1 ); clse choose y\ = y] = V(-yJj; Apply 
vj =-G(y l -y\,y 2 ) , 
v 2 = ~G(y 3 , y A ) , 

until O' i» >2> y-*' y 4* > 5 ) = O’ 1 . 0. 0, 0, y 5 l ); then go to Maneuver 3. 

Maneuver 3: Apply 

v, = - GO’i -y\,y 2 ) , 
v 2 = - G£y 3 -J 3 , y 4 ) , 

until O' 1 . > 2 . >3> > 4 . > 5 ) = 0 1 . 0. yl . 0, 0); then go to Maneuver 4. 

Maneuver 4: Apply 

Vi = - GO’i. yi) > 
v 2 = -G(y 3 ~ y 3 , y 4 ) , 

until O’ 1 . y 2 * y 3> >’4- y 5 ) = (0, 0 , > 3 , 0 , 0); then go to Maneuver 5. 

Maneuver 5: Apply 

V! = - GOi. y %) , 

v 2 = - G(y 3 ,y 4 ) , 

until Oi. 3'2» y 3- y 4» ^ 5 ) = (0- 0, 0, 0, 0); then go to Maneuver 2. 

It can be verified that the execution of Maneuver 1 transfers the initial state of the nor- 
mal form equations to the equilibrium state (0, 0, 0, 0, >>5 ), for some y 5 , in finite time. Sub- 
sequent execution of Maneuvers 2 through 5 then transfers the state (0, 0, 0, 0, y 5 ) to the ori- 
gin in finite time. This control algorithm is nonclassical and involves switching between vari- 
ous feedback functions. Justification that it stabilizes the origin of the normal form equations 
( 4 . 7 )-( 4 . 11 ) in finite time follows as a consequence of the construction procedure. Since stabil- 
ization of the normal form equations to the origin is equivalent to stabilization of the state 
equations ( 3 . 7 )-( 3 . 11 ) to the origin, we conclude that the control inputs u { and u 2 given by 
equation ( 4 . 6 ) with v t and v 2 defined by the above control algorithm stabilizes the reduced 
spacecraft dynamics described by equations (3.7)-(3. 1 1) to the equilibrium (coj, © 2 , 4*> Q. V) = 
(0, 0, 0, 0, 0) in finite time. A computer implementation of the feedback control strategy can 
be easily carried out. 
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4.2. Feedback stabilization based on rigid body rotational characteristics 


We now present an alternate discontinuous feedback control strategy for stabilizing the 
origin of equations (3.7)-(3.11) in finite time. This strategy requires that the spacecraft 
undergo a sequence of specified maneuvers and is based on the fact that rigid body rotations 
do not commute. The physical interpretation of the sequence of maneuvers that transfers any 
initial state of equation (3.7)-(3.1 1) to the origin is as follows. 






Transfer the initial state of equations (3.7)-(3.11) to any equilibrium state in finite time; 
i.e. bring the spacecraft to rest. 

Transfer the resulting state to an equilibrium state where <j> = 0 in finite time; i.e. so that 
the spacecraft is at rest with <J> = 0. 

Transfer the resulting state to an equilibrium state where <)> = 0, 0 = 0 in finite time; i.e. 
so that the spacecraft is at rest with <)) = 0, 0 = 0. 


71 


Transfer the resulting state to an equilibrium state where <j> = — , 0 = 0 in finite time; i.e. 


n 


so that the spacecraft is at rest with <J> = y, 0 = 0. 


• Transfer the resulting state to the equilibrium state (0, 0, y, 0, 0) in finite time. 

• Transfer the equilibrium state (0, 0, y, 0, 0) to the equilibrium state (0, 0, 0, 0, 0) in 
finite time. 


We now present a feedback control algorithm which stabilizes the spacecraft to the origin in 
finite time; this feedback control algorithm implements the approach just described. 


Maneuver 1. Apply 

u { = - k sign©! , 
u 2 = - k signup , 

until (©!, © 2 ) = (0. 0); then go to Maneuver 2. 
Maneuver 2: Apply 

= - G(<|>, ©!) , 

U 2 = 0 , 
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until (0)j, CO 2 , <J>) = (0, 0, 0); then go to Maneuver 3. 

Manuever 3: Apply 

«i = 0 , 

u 2 = - G(0, co 2 ) , 

until (G) lt 0 ) 2 , <{>, 0) = (0, 0, 0, 0); then go to Maneuver 4. 

Maneuver 4: Apply 

u x =- G(< j> - y, «!) , 

u 2 = 0, 

until (0)!, 0 ) 2 , 0, 0) = (0, 0, y, 0), then go to Maneuver 5. 

Maneuver 5: Apply 

u i=0, 

u 2 = - G(v, 0)2) , 

until (co lt 0 ) 2 , <|>, 0, V) = (0, 0, y, 0, 0); then go to Maneuver 6. 

Maneuver 6: Apply 

u , = - G (4>, Oj) , 

u 2 = 0 , 

until (©!, 0 ) 2 , 4>, 0, y) = (0, 0, 0, 0, 0, 0); then go to Maneuver 1. 

It can be verified that the execution of Maneuver 1 transfers the initial state of equations 
(3.7)-(3.11) to the equilibrium state (0, 0, V, 0 1 , y 1 ), for some 0 1 , 0 1 , y 1 , in finite time. Exe- 
cution of Manuever 2 then transfers the state (0, 0, <►*, 0 1 , y l ) to the state (0, 0, 0, 0 1 , y 1 ); 
execution of Manuever 3 then transfers the state (0, 0, 0, 0 1 , y 1 ) to the state (0, 0, 0, 0, y 1 ); 

execution of Manuever 4 then transfers the state (0, 0, 0, 0, y 1 ) to the state (0, 0, y , 0, y 1 ); 

execution of Manuever 5 then transfers the state (0, 0, y, 0, y 1 ) to the state (0, 0, y , 0, 0); 

finally, execution of Manuever 6 transfers the state (0, 0, y , 0, 0) to the state (0, 0, 0, 0, 0). 

This strategy is discontinuous and nonclassical in nature. A computer implementation of the 
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feedback control strategy can be easily earned out 
4.3 Comments 

We have introduced two different control laws which transfer any initial state of equa- 
tions (3.7M3.11) to the origin in finite time. Each of these control laws is in feedback form, 
since the control values depend on the current state; and each control law is discontinuous. 
The first construction procedure makes use of the nonholonomic features of the reduced 
spacecraft dynamics, while the second construction procedure uses physical insight about rigid 
body rotations. Thfe first control law constructed makes use of both control actuators simul- 
taneously, while the second control law (after Maneuver 1) uses only a single actuator at a 
time. The two discontinuous feedback control laws exhibited are illustrations of the class of 
control laws which asymptotically stabilize equations (3.7)-(3.11) to the origin. There are 
other maneuver sequences, and corresponding feedback control laws, which will also achieve 
the desired attitude stabilization of the spacecraft. But each such strategy is necessarily 
discontinuous. 

One of the advantages of the development in Sections 4.1 and 4.2 is that feedback con- 
trol strategies are constructed which guarantee attitude stabilization in a finite time. The total 
time required to complete the spacecraft reorientation is the sum of the times required to com- 
plete the sequence of maneuvers described. It should be clear that the time required to com- 
plete each maneuver depends on the single positive parameter k in the corresponding control 
law. There is a trade off between the required control levels, determined by the selection of k , 
and the resulting tim es to complete each of the maneuvers and hence the total time required to 
reorient the spacecraft In particular, the time to reorient the spacecraft from a given initial 
state to the origin can be expressed as a function of the value of the parameter k and of the 
initial state. 

We have demonstrated, by construction, the closed loop properties for the special feed- 
back control strategies presented. Our analysis was based on an ideal model assumption. 
Further robustness analysis is required to determine effects of model uncertaimties and exter- 
nal dis turbances. Unfortunately, such robustness analysis is quite difficult since the closed 
loop vector fields are necessarily discontinuous. Perhaps, feedback control strategies which 
stabilize the spacecraft attitude, different from ones presented in this paper, would provide 
improved closed loop robustness. These issues are to be studied in future research. 
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5. Simulation 

We illustrate the results of the paper using an example. Consider a rigid spacecraft with 
no control torque about the third principal axis and two control torques, generated by momen- 
tum wheel actuators, are applied about the other two principal axes. Therefore the vectors b j 
and bi are given by b\ = (1, 0, 0) T , b 2 = (0, 1, 0) T . For our simulation, we use the space- 
craft parameters used in Ref. 2. The mass of the spacecraft, m j, is 500 Kg, and the masses of 
the momentum wheels, m 2 and m3, are each 5 Kg. The center of mass of the momentum 
wheels are located at a distance 0.2 m from the center of mass of the spacecraft, i.e., d x - d 2 
= 0.2 m. The moment of inertia of the wheels about its axis of rotation is 0.5 Kg.m 2 , i.e., 
j j = j 2 = 0.5. The inertia tensor of the spacecraft and the two momentum wheels are 

1 1 = diag (86.215, 85.07, 113.565) Kg.m 2 , 

1 2 = diag (0.5, 0.25, 0.25) Kg.m 2 , 

1 2 = diag (0.25, 0.5, 0.25) Kg.m 2 . 

Using these parameters, the inertia matrix J can be calculated which equals 
J = diag(86.7, 85.5, 1 14.5) Kg.m 2 , 

approximately. The complete dynamics of the spacecraft system defined by equations (3.1)- 
(3.4) is not controllable, but we consider the restriction that the angular momentum vector 
H = 0. Consequently, we are interested in stabilizing the reduced spacecraft dynamics 
described by equations (3.7)-(3.11) to the equilibrium (coj, o>2> <i>. 6, V) = (0, 0, 0, 0, 0). 
The spacecraft is initially at rest (i.e., co® = 00° = 0) with an initial orientation given by the 
Euler angles <)> 0 = re, 0° = 0.25jt and \jr° = - 0.57t. 

First, a computer implementation of the feedback control algorithm specified in Section 
4.1 was used to stabilize the spacecraft to the origin. The value of the gain k was chosen as 
1. The time responses of the Euler angles, angular velocities and the control torques are 
shown in Fig. 1, Fig. 2 and Fig. 3 respectively. After a total maneuver time of 11.77 seconds, 
( o 1 = c0 2 = 4> = e = y = O. Next, a computer implementation of the feedback control algo- 
rithm specified in Section 4.2 was used to stabilize the spacecraft to the origin. The value of 
the gain k was chosen as 1. The time responses of the Euler angles, angular velocities and 
the control torques are shown in Fig. 4, Fig. 5 and Fig. 6 respectively. After a total maneuver 
time of 13 seconds, co 1 = o)2 = 0 = 0 = y = O. 
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6. Conclusion 

The attitude stabilization problem of a spacecraft using control torques supplied by two 
momentum wheel actuators about axes spanning a two dimensional plane orthogonal to a prin- 
cipal axis has been considered. The complete spacecraft dynamics are not controllable. How- 
ever, the spacecraft dynamics are small time locally controllable in a reduced sense. The 
reduced spacecraft dynamics cannot be asymptotically stabilized using time-invariant continu- 
ous feedback, but discontinuous feedback control strategies have been constructed which sta- 
bilizes the spacecraft (in the reduced sense) to an equilibrium attitude in finite time. The 
results of the paper show that although classical nonlinear control techniques do not apply, it 
is possible to construct control laws based on the particular spacecraft dynamics. 
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Fig. 1: Euler Angles 



Fig. 2: Anglular Velocities 






Fig. 3: Control Torques 



Fig. 4: Euler Angles 
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Fig. 5: Anglular Velocities 



Fig. 6: Control Torques 







